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ABSTRACT 

The primary nit rahigh energy particles which produce giant extensive air 
showers in the Earth’s atmosphere present an intriguing mystery from two points 
of view: (1) How are these particles produced with such astounding energies, eight 
orders of magnitude higher than those produced by the best man-made terrestrial 
accelerators? (2) Since they are most likely extragalactic in origin, how do they 
reach us from extragalactic distances without suffering the severe losses expected 
from interactions with the 2.7 K thermal cosmic background photons - the so- 
called GZI< effect? 

The answers to these questions may involve new physics: violations of spe- 
cial relativity, grand unification theories, and quantum gravity theories involving 
large extra dimensions. They may involve new astrophysical sources, ”zeva- 
trons”. Or some heretofore totally unknown physics or astrophysics may hold 
the answer. I will discuss here the mysteries involving the production and extra- 
galactic propagation of ultrahigh energy cosmic rays and some suggested possible 
solutions. 


Subject headings: ultrahigh energy cosmic rays, active galactic nuclei, gamma-ray 
bursts, topological defects, grand unification 


1. Introduction 

About once per century per km 2 of the Earth’s surface, a giant shower of charged 
particles produced by a primary particle with an energy greater than or equal to 16 joules 
(100 EeV = lO 20 eV) plows through the Earth’s atmosphere. The showers which they 
produce can be detected by arrays of scintillators on the ground; they also announce their 
presence by producing a trail of ultraviolet flourescent light, exciting the nitrogen atoms 
in the atmosphere. The existence of such showers has been known for almost four decades 
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(Linsley 1963). The number of giant air showers detected from primaries of energy greater 
than 100 EeV has grown into the double digits and may grow into the hundreds as new 
detectors such as the “Auger” array and the “EUSO” (Extreme Universe Space Observatory) 
and “OWL” (Orbiting Wide-Angle Light Collectors) satellite detectors come on line. These 
phenomena present an intriguing mystery from two points of view: (1) How are particles 
produced with such astounding energies, eight orders of magnitude higher than are produced 
by the best man-made terrestrial accelerators? (2) Since they are most likely extragalactic 
in origin, how do they reach us from extragalactic distances without exhibiting the predicted 
cutoff from interactions with the 2.7I< cosmic background radiation? In these lectures, I will 
consider possible solutions to this double mystery. 


2. The Data 

Figure 1 shows the published data (as of this writing) on the ultrahigh energy cosmic 
ray spectrum from the Fly’s Eye and AG AS A detectors. 1 Other data from Havera Park and 
Yakutsk may be found in the review' by Nagano and Watson (2000) are consistent with 
Figure 1. Additional data are now being obtained by the HiRes detector array and should 
be available in the near future (T. Abu-Zayyad, et al. , in preparation). 

For air showers produced by primaries of energies in the 1 to 3 EeV range, Hayashida, 
et al. (1999) have found a marked directional anisotropy with a 4.5<r excess from the galactic 
center region, a 3.9cr excess from the Cygnus region of the galaxy, and a 4.0<r deficit from 
the galactic anticenter region. This is strong evidence that EeV cosmic rays are of galactic 
origin. A galactic plane enhancement in EeV events was also reported by the Fly’s Eye 
group (Dai, et al. 1999). 

As shown in Figure 2, at EeV energies, the primary particles appear to have a mixed or 
heavy composition, trending toward a light composition in the higher energy range around 
30 EeV (Bird, et al. 1993; Abu-Zayyad, et al. 2000). This trend, together with evidence of 
a flattening in the cosmic ray spectrum on the 3 to 10 EeV energy range (Bird, et al. 1994; 
Takeda et al. 1998) is evidence for a new component of cosmic rays dominating above 10 
EeV energy. 

The apparent isotropy (no galactic-plane enhancement) of cosmic rays above 10 EeV 
{e.g. Takeda, et al. 1999), together with the difficulty of confining protons in the galaxy at 


'The AG ASA data have been reanalysed and the number of events determined to be above 100 EeV has 
been lowered to eight. (Teshima, private communication.) 
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10 to 30 EeV energies, provide significant reasons to believe that the cosmic-ray component 
above 10 EeV is extragalactic in origin. As can be seen from Figure 1, this extragalactic 
component appears to extend to an energy of 300 EeV. Extention of this spectrum to higher 
energies is conceivable because such cosmic rays, if they exist, would be too rare to have 
been seen with present detectors. We will see in the next section that the existence of 300 
EeV cosmic rays gives us a new mystery to solve. 


3. The GZK Effect 

Thirty seven years ago, Penzias and Wilson (1965) reported the discovery of the cosmic 
3K thermal blackbody radiation which was produced very early on in the history of the 
universe and which led to t he undisputed acceptance of the “big bang” theory of the origin 
of the universe. Much more recently, the Cosmic Background Explorer (COBE) satellite 
confirmed this discovery, showing that the cosmic background radiation (CBR) has the 
spectrum of the most perfect thermal blackbody known to man. COBE data also showed 
that this radiation (on angular scales > 7°) was isotropic to a part in 10 5 (Mather et al. 
1994). The perfect thermal character and smoothness of the CBR proved conclusively that 
this radiation is indeed cosmological and that, at the present time, it fills the entire universe 
with a 2.7K spectrum of radio to far-infrared photons with a density of ~ 400 cm 3 . 

Shortly after the disc overy of the CBR, Greisen (1966) and Zatsepin and Kuz'min 
(1966) predicted that pion-producing interactions of ultrahigh energy cosmic ray protons 
with CBR photons of target density ~ 400 cm -3 should produce a cutoff in their spectrum 
at energies greater than ~ 50 EeV. This predicted effect has since become known as the GZK 
(Greisen-Zatsepin-Kuz'min) effect, hollowing the GZK papers, Stecker (1968) utilized data 
on the energy dependence of the photomeson production cross sections and inelasticities to 
calculate the mean energy loss time for protons propagating through the CBR in intergalactic 
space as a function of energy. Based on his results, Stecker (1968) then suggested that the 
particles of energy above the GZK cutoff energy (hereafter referred to as trans-GZK particles) 
must be coming from within the “Local Supercluster” of which we are a part and which is 
centered on the Virgo Cluster of galaxies. Thus, the “GZK cutoff” is not a true cutoff, but 
a supression of the ultrahigh energy cosmic ray flux owing to a limitation of the propagation 
distance to a few tens of Mpc. 

The actual position of the GZK cutoff can differ from the 50 EeV predicted by Greisen. 
In fact, there could actually be an enhancement at or near this energy owing to a “pileup” 
of cosmic rays starting oui at higher energies and crowding up in energy space at or below 
the predicted cutoff energy (Puget, Stecker and Bredkamp 1976; Hill and Schramm 1985; 
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Fig. 1. — The ultrahigh energy cosmic ray spectrum data from Fly’s Eye and AG AS A. 



Fig. 2. — Average depth of shower maximum (A max ) vs. energy compared to the calculated 
values for protons (upper curves) and Fe primaries (lower curves) (from Gaisser 2000; see 
references therein). 
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Berezinsky and Grigor’eva 1988; Stecker 1989; Stecker and Salamon 1999). The existence 
and intensity of this predicted pileup depends critially on the flatness and extent of the 
source spectrum, (z. 6., the number of cosmic rays starting out at higher energies), but if its 
existence is confirmed in the future by more sensitive detectors, it would be evidence for the 
GZK effect. 

Scully and Stecker (2002) have determined the GZK energy, defined as the energy for a 
flux decrease of 1/e, as a i unction of redshift . At high redshifts, the target photon densitv 
increases by (1 + zf and both the photon and initial cosmic ray energies increase by (1 + z). 
The results obtained by Scully and Stecker are shown in Figure 3. 


4. Acceleration and Zevatrons: The “Bottom Up” Scenario 

The apparent lack of a GZK cutoff has led theorists to go on a hunt for neat by ze- 
vatrons”, z. e,, astrophysics! sources which can accelerate particles to energies 0(1 ZeV 
10 21 eV). 

In most theoretical work in cosmic ray astrophysics, it is generally assumed that the 
diffusive shock acceleration process is the most likely mechanism for accelerating particles 
to high energy. (See, e.$., Jones (2000) and references therein.) In this case, the maximum 
obtainable energy is given by E max = keZ{u/c)BL , where u < c is the shock speed, eZ is 
the charge of the particle being accelerated, B is the magnetic field strength, L is the size of 
the accelerating region and the numerical parameter k = 0(1) ( Drury 1994). Taking k = 1 
and u = c, one finds 

Emax = 0.9 Z(BL) 

with E in EeV, B in //G and R in kpc. This assumes that particles can be accelerated 
efficiently up until the moment when they can no longer be contained by the source, i. t. 
until their gyroradius becomes larger than the size of the source. Hillas (1984) used this 
relation to construct a plot, of B vs. L for various candidate astrophysical objects. A Hillas 
plot” of this kind, recently constructed by Olinto (2000), is shown in Figure 4. 

Given the relationship between E max and BL as shown in Figure 4, there are not too 
many astrophysical candidates for zevatrons. Of these, galactic sources such as white dwarfs, 
neutron stars, pulsars, and magnetars can be ruled out because their galactic distribution 
would lead to anisotropies above 10 EeV which would be similar to those observed at lower 
energies by Hayashida et al (1999), and this is not the case. Perhaps the most promising 



- 6 - 



t + z 


Fig. 3. — The GZK cutoff energy versus redshift (Scully and Stecker 2002). 



Fig. 4. — A “Hillas Plot” showing potential astrophysical zevatrons (from Olinto 2000). The 
lines are for B vs. L for E max = 0.1 ZeV for protons and iron nuclei as indicated. 
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potential zevatrons are radio lobes of strong radio galaxies (Biermann and Strittmatter 
(1987) . The trick is that such sources need to be found close enough to avoid the GZK 
cutoff ( e.g ., Elbert and Sommers 1995). Biermann has further suggested that the nearby 
radio galaxy M87 may be t he source of the observed trans-GZK cosmic rays (see also Stecker 
1968; Farrar and Piran 2000). Such an explanation would require one to invoke magnetic 
field configurations capable of producing a quasi-isotropic distribution of > 10 2 eV protons, 
making this hypothesis questionable. However, if the primary particles are nuclei, it is easier 
to explain a radio galaxy origin for the two highest energy events (Stecker and Salamon 1999; 
see section 6). 

It has also been suggested that since all large galaxies are suspected to harbor super- 
massive black holes in their centers which may have once been quasars, fed by accretion 
disks which are now used up, that nearby quasar remnants may be the searched-for zeva- 
trons (Boldt and Ghosh 1999; Boldt and Lowenstein 2000) . This scenario also has potential 
theoretical problems and needs to be explored further. In particular, it has been shown 
that black holes which an' not accreting plasma cannot possess a large scale magnetic field 
with which to accelerate particles to relativistic energies (Ginzburg and Ozernoi 1964; Krolik 
1999; Jones 2000). Observational evidence also indicates that the cores of weakly active 
galaxies have low magnetic fields (Falcke 2001 and references therein). Another proposed 
zevatron, the 7-ray burst, is discussed in the next section. 


5. Gamma-Ray Burst Zevatrons and the GZK Problem 

In 1995, it was suggested that cosmological 7-ray bursts (GRBs) were the source of the 
highest energy cosmic rays (Waxman 1995; Vietri 1995). It was suggested that if these 
objects emitted the same amount of energy in ultrahigh energy (~ 10 14 MeV) cosmic rays as 
in ~ MeV photons, there would be enough energy input of these particles into intergalactic 
space to account for the 3 observed flux. At that time, it was assumed that the GRBs were 
distributed uniformly, independent of redshift. 

In recent years, X-rav, optical, and radio afterglows of about a dozen GRBs have been 
detected leading to the subsequent identification of the host galaxies of these objects and 
consequently, their redshifts. The host galaxies of GRBs appear to be sites of active star 
formation. The colors and morphological types of the host galaxies are indicative of ongoing 
star formation, as is the detection of Lya and [Oil] in several of these galaxies. Further 
evidence suggests that bursts themselves are directly associated with star forming regions 
within their host galaxies their positions correspond to regions having significant hydrogen 
column densities with evidence of dust extinction. It now seems more reasonable to assume 
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that a more appropriate redshift distribution to take for GRBs is that of the average star 
formation rate. 

To date, some 14 GRBs afterglows have been detected with a subsequent identification 
of their host galaxies. As of this writing, 13 of the 14 are at moderate to high redshifts with 
the highest one (GRB000131) lying at a redshift of 4.50 (Andersen, et ai 2000). 

A good argument in favor of strong redshift evolution for the frequency of occurrence of 
the higher luminosity GRBs has been made by Mao and Mo (1998), based on the nature of 
the host galaxies. Other recent analyses have also favored a GRB redshift distribution which 
follows the strong redshift evolution of the star formation rate (Schmidt 1999; Fenimore 
and Ramirez-Ruiz 2000). If we thus assume a redshift distribution for the GRBs which 
follows the star formation rate, being significantly higher at higher redshifts, GRBs fail by at 
least an order of magnitude to account for the observed cosmic rays above 100 EeV (Stecker 
2000). If one wishes to account for the GRBs above 10 EeV, this hypothesis fails by two to 
three orders of magnitude (Scully and Stecker 2002). Even these numbers are most likely 
too optimistic, since they are based on the questionable assumption of the same amount of 
GRB energy being put into ultrahigh energy cosmic rays as in ~ MeV photons. 

Figure 5, from Scully and Stecker, (2002) shows the form of the cosmic ray spectrum to 
be expected from sources with a uniform redshift distribution and sources which follow the 
star formation rate. The required normalization and spectral index determine the energy 
requirements of any cosmological sources which are invoked to explain the observations. 
Pileup effects and GZK cutoffs are evident in the theoretical curves in this figure. As can 
be seen in Figure 5, the present data appear to be statistically consistent with either the 
presence or absence of a pileup effect. Future data with much better statistics are required 
to determine such a spectral structure. 

An unusual nearby Type Ic supernova, SN 1998bw, has been identified as the nearby 
source of a low luminosity burst, GRB980425, with an energy release which is orders of 
magnitude smaller than that for a typical cosmological GRB. Norris (2002) has given an 
analysis of the luminosities and space densities of such nearby low luminosity long-lag GRB 
sources which are identified with Type I supernovae. For these sources, he finds a rate per 
unit volume of 7.8 x 10” 7 Mpc“ 3 yr _1 and an average (isotropic) energy release per burst of 
1.3 x 10 49 erg over the energy range from 10 to 1000 keV. The energy release per unit volume 
is then ~ 10 43 erg Mpc^yr -1 . This rate is more than an order of magnitude below the rate 
needed to account for the cosmic rays with energies above 10 EeV. 
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Fig. 5. — Predicted spectra for cosmic ray protons as compared with the data. The middle 
curve and lowest curve assume an £~ 2 ' 75 source spectrum with a uniform source distribution 
and one that follows the : distribution of the star formation rate respectively. The upper 
curve is for an E ~ 235 source spectrum which requires an order of magnitude more energy 
input and exhibits a “pileup effect” (Scully and Stecker 2002). 
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6. The Heavy Nuclei Origin Scenario 

A more conservative hypothesis for explaining the trans-GZK events is that they were 
produced by heavy nuclei. Stecker and Salamon (1999) have shown that the energy loss time 
for nuclei starting out as Fe is longer than that for protons for energies up to a total energy 
of ~300 EeV (see Figure 6). 

Stanev, et al. (1995) and Biermann (1998)1 have examined the arrival directions of the 
highest energy events. They point out that the ~ 200 EeV event is within 10° of the direction 
of the strong radio galaxy NGC 315. This galaxy lies at a distance of only ~ 60 Mpc from 
us. For that distance, the results of Stecker and Salamon (1999) indicate that heavy nuclei 
would have a cutoff energy of ~ 130 EeV, which may be within the uncertainty in the energy 
determination for this event. The ~300 EeV event is within 12° of the direction of the strong 
radio galaxy 3C134. The distance to 3C134 is unfortunately unknown because its location 
behind a dense molecular cloud in our own galaxy obscures the spectral lines required for 
a measurement of its redshift. It may be possible that either cosmic ray protons or heavy 
nuclei originated in these sources and produced the highest energy air shower events. 

An interesting new clue that we may indeed be seeing heavier nuclei above the proton- 
GZI< cutoff comes from a very recent analysis of inclined air showers above 10 EeV energy 
(Ave, et al. 2000). These new results favor proton primaries below the p-GZK cutoff energy 
but they appear to favor a heavier composition above the p-GZK cutoff energy. It will be 
interesting to see what future data from much more sensitive detectors will tell us. 


7. Top-Down Scenarios: “FYaggers” 

A way to avoid the problems with finding plausible astrophysical zevatrons is to start 
at the top, he., the energy scale associated with grand unification, supersymmetric grand 
unification or its string theory equivalent. 

The modern scenario for the early history of the big bang takes account of the work of 
particle theorists to unify the forces of nature in the framework of Grand Unified Theories 
(GUTs) (e.g., Georgi and Glashow 1974). This concept extends the very successful work of 
Nobel Laureates Glashow, Weinberg, and Salam in unifying the electromagnetic and weak 
nuclear forces of nature (Glashow 1960; Weinberg 1967; Salam 1968). As a consequence 
of this theory, the electromagnetic and weak forces would have been unified at a higher 
temperature phase in the early history of the universe and then would have been broken 
into separate forces through the mechanism of spontaneous symmetry breaking caused by 
vacuum fields which are known as Higgs fields. 
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In GUTs, this same paradigm is used to infer that the electroweak force becomes unified 
with the strong nuclear force at very high energies of ~ 10 24 eV which occurred only ~ 
10 -35 seconds after the big bang. The forces then became separated owing to interactions 
with the much heavier mass scale Higgs fields whose symmetry was broken spontaneously. 
The supersymmetric GUTs (or SUSY GUTs) provide an explanation for the vast difference 
between the two unification scales (known as the “Hierarchy Problem”) and predict that the 
running coupling constants which describe the strength of the various forces become equal 
at the SUSY GUT scale of ~ 10 24 eV (Dimopoulos, Raby and Wilczek 1982). 


7.1. Topological Defects 

The fossil remnants of this unification are predicted to be very heavy topological defects 
in the vacuum of space caused by misalignments of the heavy Higgs fields in regions which 
were causally disconnected in the early history of the universe. These are localized regions 
where extremely high densities of mass-energy are trapped. Such defects go by designations 
such as cosmic strings, monopoles, walls, necklaces (strings bounded by monopoles), and 
textures, depending on their geometrical and topological properties. Inside a topological 
defect, the vestiges of the early universe may be preserved to the present day. The general 
scenario for creating topological defects in the early universe was suggested by Kibble (1976). 

Superheavy particles or topological structures arising at the GUT energy scale M > 10 23 
eV can decay or annihilate to produce “X-particles” (GUT scale Higgs particles, superheavy 
fermions, or leptoquark bosons of mass M.) In the case of strings this could involve mecha- 
nisms such as intersecting and intercommuting string segments and cusp evaporation. These 
X-particles will decay to produce QCD fragmentation jets at ultrahigh energies, so I will 
refer to them as “fraggers”. QCD fraggers produce mainly pions, with a 3 to 10 per cent 
admixture of baryons, so that generally one can expect them to produce at least an order 
of magnitude more ultrahigh energy 7-rays and neutrinos than protons. The same general 
scenario would hold for the decay of long-lived superheavy dark matter particles (see section 
7.3), which would also be fraggers. It has also been suggested that the decay of ultraheavy 
particles from topological defects produced in SUSY-GUT models which can have an addi- 
tional soft symmetry breaking scale at TeV energies (“flat SUSY theories”) may help explain 
the observed 7-ray background flux at energies ~ 0.1 TeV (Bhattacharjee, Shafi and Stecker 
1998). 

The number of variations and models for explaining the ultrahigh energy cosmic rays 
based on the GUT or SUSY GUT scheme (which have come to be called “top-down” models) 
has grown to be enormous and I will not attempt to list all of the numerous citations involved. 
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Fortunately, Bhattacharjee and Sigl (2000) have recently published an extensive review with 
over 500 citations and I refer the reader to this review for further details of “top-down” models 
and references. The important thing to note here is that, if the implications of such models 
are borne out by future cosmic ray data, they may provide our first real evidence for GUTs. 


7.2. “Z-bursts” 

It has been suggested that ultra-ultrahigh energy (9(10 ZeV) neutrinos can produce 
ultrahigh energy Z 0 fraggrrs by interactions with 1.9K thermal CBR neutrinos (Weiler 1982, 
Fargion, et al 1999; Weiler 1999), resulting in “Z-burst” fragmentation jets, again producing 
mostly pions. This will oc< ur at the resonance energy E res = 4[m l/ (eV)]‘ 1 ZeV . A typical Z 
boson will decay to produce ^2 nucleons, ~20 7-rays and ~ 50 neutrinos, 2/3 of which are 

*V s - 

If the nucleons which are produced from Z-bursts originate within a few tens of Mpc 
of us they can reach us, even though the original ~ 10 ZeV neutrinos could have come 
from a much further distance. It has been suggested that this effect can be amplified if our 
galaxy has a halo of neutrinos with a mass of tens of eV (Fargion, Mele and Salis 1999; Weiler 
1999). However, a neutrino mass large enough to be confined to a galaxy size neutrino halo 
(Tremaine and Gunn 1979) would imply a hot dark matter cosmology which is inconsistent 
with simulations of galaxy formation and clustering (e.#., Ma and Bertschinger 1994) and 
with angular fluctuations in the CBR. (Another problem with halo fraggers is discussed 
below in section 7.4) A mixed dark matter model with a lighter neutrino mass (Shaft and 
Stecker 1984) produces predicted CBR angular fluctuations (Schaefer, Shaft and Stecker 
1989) which are consistent with the Cosmic Background Explorer data (Wright 1992). In 
such a model, neutrinos would have density fluctuations on the scale of superclusters, which 
would still allow for some amplification (Weiler 1999) . The tritium decay spectral endpoint 
limits on the mass of the electron neutrino (Weinheimer, et al. 1999), together with the 
very small neutrino flavor mass differences indicated by the atmospheric and solar neutrino 
oscillation results (Ahmad, et al. 2002) constrains all neutrino flavors to have masses in the 
range 0(eW) or less. This is much too small a mass for neutrinos to to be confined to halos 
of individual galaxies. 

The basic general problem with the Z-burst explanation for the trans-GZK events is 
that one needs to produce 10 ZeV neutrinos. If these are secondaries from pion production, 
this implies that the primary protons which produce them must have energies of hundreds of 
ZeV! Since we know of no astrophysical source which would have the potential of accelerating 
particles to energies even an order of magnitude lower (see section 4), a much more likely 
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scenario for producing 10 ZeV neutrinos would be by a top-down process. The production 
rate of neutrinos from such processes is constrained by the fact that the related energy 
release into electromagnetic cascades which produce GeV range 7-rays is limited by the 
satellite observations (see the review by Bhattacharjee and Sigl 2000). This constraint, 
together with the low probability for Z-burst production, relegates the Z-burst phenomenon 
to a minor secondary role at best. 


7.3. Ultraheavy Dark Matter Particles: “Wimpzillas” 

The idea has been suggested that the dark matter which makes up most of the gravi- 
tating mass in the universe could consist of ultraheavy particles produced by non-thermal 
processes in the early big-bang (Berezinsky et ai 1997; Kuz’min and Rubakov 1998; Blasi 
et ai 2002). (See also the paper of Rocky Kolb in these proceedings.) The annihilation 
or decay of such particles in a dark matter halo of our galaxy would then produce ultra- 
high energy nucleons which would not be attenuated at trans-GZK energies owing to their 
proximity. 


7.4. Halo Fraggers and the Missing Photon Problem 

Halo fragger models such as Z-burst and ultraheavy halo dark matter (“wimpzilla”) 
decay or annihilation, as we have seen, will produce many more photons than protons. These 
ultrahigh energy photons can reach the Earth from anywhere in a dark matter galactic halo, 
because, as shown in Figure 7, there is a “mini-window” for the transmission of ultrahigh 
energy cosmic rays between ~ 0.1 and ~ 10 6 EeV. 

Photon-induced giant air showers have an evolution profile which is significantly different 
from nucleon-induced showers because of the Landau-Pomeranchuk-Migdal (LPM) effect 
(Landau and Pomeranchuk 1953; Migdal 1956) and because of cascading in the Earth’s 
magnetic field (Cillis, et ai 1999) (see Figure 7). By taking this into account, Shinozaki, 
et ai (2002) have used the AGASA data to place upper limits on the photon composition 
of their UHECR showers. They find a photon content upper limit of 28% for events above 
10 EeV and 67% for events above 30 EeV at a 95% confidence level with no indication of 
photonic showers above 100 EeV. A recent reanalysis of the ultrahigh energy events observed 
at Haverah Park by Ave, et ai (2002) indicates that less than half of the events (at 95% 
confidence level) observed above 10 and 40 EeV are 7-ray initiated. An analysis of the 
highest energy Fly’s Eye event (E= 300 EeV) (Halzen and Hooper 2002) shows it not to be 
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of photonic origin, as indicated in Figure 8. In addition, Shinozaki, et al. (2002) have found 
no indication of departures from isotropy as would be expected from halo fragger photonic 
showers, this admittedly with only 10 events in their sample. 


8. Other New Physics Possibilities 

The GZK cutoff problem has stimulated theorists to look for possible solutions involving 
new physics. Some of these involve (A) a large increase in the neutrino-nucleon cross section 
at ultrahigh energies, (B) new particles, and (C) a small violation of Lorentz Invariance (LI). 


8.1. Increasing the Neutrino-Nucleon Cross Section at Ultrahigh Energies 

Since neutrinos can travel through the universe without interacting with the 2./ K CBR, 
it has been suggested that if the neutrino-nucleon cross section were to increase to hadronic 
values at ultrahigh energies, they could produce the giant air showers and account for the 
observations of showers above the proton-GZK cutoff. Several suggestions have been made 
for processes that can enhance the neutrino-nucleon cross section at ultrahigh energies. These 
suggestions include composite models of neutrinos (Domokos and Nussinov 1987; Domokos 
and Kovesi-Domokos 1988) , scalar leptoquark resonance channels (Robinett 1988) and the 
exchange of dual gluons (Bordes, et al. 1998). Burdman, Halzen and Ghandi (1998) have 
ruled out a fairly general class of these types of models, including those listed above, by 
pointing out that in order to increase the neutrino-nucleon cross section to hadronic values 
at ~ lO 20 eV without violating unitarity bounds, the relevant scale of compositeness or 
particle exchange would have to be of the order of a GeV, and that such a scale is ruled out 
by accelerator experiments. 

More recently, the prospect of enhanced neutrino cross sections has been explored in 
the context of extra dimension models. Such models have been suggested by theorists to 
unify the forces of physics since the days of Kaluza (1921) and Klein (1926). In recent 
years, they have been invoked by string theorists and by other theorists as a possible way 
for accounting for the extraordinary weakness of the gravitational force, or, in other words, 
the extreme size of the Pianck mass (Arkani-Hamed, Dimopoulos and Dvali 1999; Randall 
and Sundrum 1999). These models allow the virtual exchange of gravitons propagating 
in the bulk (i.e. in the space of full extra dimensions) while restricting the propagation of 
other particles to the familiar four dimensional space-time manifold. It has been suggested 
that in such models, a{u M) ~ [£„/(! 0 2O eV)] mb (Nussinov and Schrock 1999; Jain, et al. 
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Fig. 7. — The mean free path for ultrahigh energy 7-ray attenuation vs. energy. The curve 
for electron-positron pair production off the cosmic background radiation (CBR) is based 
on Gould and Schreder (1966). The two estimates for pair production off the extragalactic 
radio background are from Protheroe and Biermann (1996). The curve for double pair 
production is based on Brown, et al. (1973). The physics of pair production by single 
photons in magnetic fields is discussed by Erber (1966). This process eliminates all photons 
above ~ 10 24 eV and produces a terrestrial anisotropy in the distribution of photon arrival 
directions above ~ 10 19 eV. 
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2000); see also Domokos and Kovesi-Domokos 1999). It should be noted that a cross 
section of ~ 100 mb would be necessary to approach obtaining consistency with the air 
shower profile data. Other scenarios involve the neutrino-initiated atmopheric production 
of black holes (Anchordoqui, et al. 2002) and even higher dimensional extended objects, 
p-dimensional branes called u p-branes” (Ahn, Cavalgia and Olinto 2002; Anchordoqui, Feng 
and Goldberg 2002). Such interactions, in principle, can increase the neutrino total 
atmospheric interaction cross section by orders of magnitude above the standard model 
value. However, as discussed by Anchordoqui, Feng and Goldberg 2002, , sub-mm gravity 
experiments and astrophysical constraints rule out total neutrino interaction cross sections 
as large as 100 mb as would be needed to fit the trans-GZK energy air shower profile data. 
Nonetheless, extra dimension models still may lead to significant increases in the neutrino 
cross section, resulting in moderately penetrating air showers. Such neutrino-induced showers 
should also be present at somewhat lower energies and provide an observational test for extra 
dimension TeV scale gravity models (Anchordoqui, et ai 2001; Tyler, Olinto and Sigl 2001). 
As of this writing, no such showers have been observed, putting an indirect constraint on 
fragger scenarios with TeY gravity models. 


8.2. New Particles 

The suggestion has also been made that new neutral particles containing a light gluino 
could be producing the trans-GZK events (Farrar 1996; Cheung, Farrar and Kolb 1998). 
While the invocation of such new particles is an intriguing idea, it seems unlikely that 
such particles of a few proton masses would be produced in copious enough quantities in 
astrophysical objects without being detected in terrestrial accelerators. Also there are now 
strong constraints on gluinos ( Alavi-Harati, et al. 1999). One should note that while it is 
true that the GZK threshold for such particles would be higher than that for protons, such 
is also the case for the more prosaic heavy nuclei (see section 7). In addition, such neutral 
particles cannot be accelerated directly, but must be produced as secondary particles, making 
the energetics reqirements more difficult. 


8.3. Breaking Lorentz Invariance 

With the idea of spontaneous symmetry breaking in particle physics came the suggestion 
that Lorentz invariance (I I) might be weakly broken at high energies (Sato and Tati 1972). 
Although no real quantum theory of gravity exists, it was suggested that LI might be broken 
as a consequence of such a theory ( Amelino-Camilia et al. 1998). A simpler formulation 
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for breaking LI by a small first order perturbation in the electromagnetic Lagrangian which 
leads to a renormalizable treatment has been given by Coleman and Glashow (1999). Using 
this formalism, these authors have shown than only a very tiny amount of LI symmetry 
breaking is required to avoid the GZK effect by supressing photomeson interactions between 
ultrahigh energy protons and the CBR. This LI breaking amounts to a difference of O(10 -23 ) 
between the maximum proton pion velocities. By comparison, Stecker and Glashow (2001) 
have placed an upper limit of 0( 10 -13 ) on the difference between the velocities of the electron 
and photon, ten orders of magnitude higher than required to eliminate the GZK effect. 


9. Is the GZK Effect All There Is? 

There is a remaining “dull” possibility. Perhaps the GZK effect is consistent with the 
data and is all there is at ultrahigh energies. The strongest case for trans-GZK physics 
comes from the AGASA results. The AGASA group, which reported up to 17 events with 
energy greater than or equal to ~ 100 EeV (Sasaki, et al. 2001), has now lowered this 
number to 8 (see footnote 1). However, the HiRes Group have not confirmed the AGASA 
results, implying lower fluxes of cosmic rays above ~ 100 EeV (T. Abu-Zayyad, et al. , in 
preparation; P. Sokolsky and E.C. Loh, private communication). Even if the GZK effect is 
seen, top-down scenarios predict the reemergence of a new component at even higher energies 
(Aharonian, Bhattacharjee and Schramm 1992; Bhattacharjee and Sigl 2000). 

The AGASA data indicate a significant deviation from pure GZK even if the source 
number is weighted like the local galaxy distribution (Blanton, et al. 2001) In addition 
to this discrepency, the fact that a flourescence detector, Fly’s Eye, reported the highest 
energy event yet seen, it viz., E ~ 300 EeV, makes the experimental situation interesting 
enough to justify both more sensitive future detectors and the exploration of new physics 
and astrophysics. 


10. Signatures 

Future data which will be obtained with new detector arrays and satellites (see next sec- 
tion) will give us more clues relating to the origin of the trans-GZK events by distinguishing 
between the various hypotheses which have been proposed. 

A zevatron origin (“bottom-up” scenario) will produce air-showers primarily from pri- 
maries w'hich are protons or heavier nuclei, with a much smaller number of neutrino-induced 
showers. The neutrinos will be secondaries from the photomeson interactions which produce 
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the GZI\ effect (Stecker 1973; 1979; Engel, Seckel and Stanev 2001 and references therein). 
In addition, zevatron events may cluster near the direction of the sources. 

A “top-down” (GUT) origin mechanism will not produce any heavier nuclei and will 
produce at least an order of magnitude more ultrahigh energy neutrinos than protons. 
Thus, it will be important to look for the neutrino-induced air showers which are expected 
to originate much more deeply in the atmosphere than proton-induced air showers and are 
therefore expected to be mostly horizontal showers. Looking for these events can most easily 
be done with a satellite array which scans the atmosphere from above (see next section). As 
we discussed, the “top dov\ n” model also produces a large ratio of ultrahigh energy photons to 
protons. This was suggested as a signature of top-down models by Aharonian, Bhatacharjee 
and Schramm (1992). However, the mean free path of these photons against pair-production 
interactions with extragalactic low frequency radio photons from radio galaxies is only a few 
Mpc at most (Protheroe and Biermann 1996). The subsequent electromagneticcasca.de and 
synchrotoron emission of the high energy electrons produced in the cascade dumps the energy 
of these particles into much lower energy photons (Wdowczyk, Tkaczyk and Wolfendale 1972; 
Stecker 1973) . The photon-proton ratio, however, is an effective tool for testing halo fragger 
models (see section 7.4). 

Another characteristic which can be used to distinguish between the bottom-up and 
top-down models is that (he latter will produce much harder spectra. If differential cosmic 
ray spectra are parametrized to be of the form F oc E ~^ , then for top-down models T < 2, 
whereas for bottom-up models T > 2. Also, because of the hard source spectrum in the 
“top-down” models, they should exhibit both a GZK suppression and a pileup just before 
the GZK energy. 

If Lorentz invariance breaking is the explanation for the missing GZK effect, the actual 
absence of photomeson interactions should result the absence of a pileup effect as well. 


11. Present and Future Detectors 

Of the ground-based ultra-high energy arrays, the AGASA array of particle detectors 
in Japan is continuing to obtain data on ultrahigh energy cosmic ray-induced air showers. 
Its aperture is 200 km 2 sr. 

The HiRes array is operating and will soon be publishing data. This array is an extension 
of the Fly’s Eye which pioneered the technique of measuring the atmospheric fluorescence 


2 For a discussion of ultrahigh energy neutrino astrophysics, see Cline and Stecker (2000). 
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Fig. 8. — The composite atmospheric shower profile of a 300 EeV photon-induced shower 
calculated with the Bethe-Heitler (solid) electromagnetic cross section and with the LPM 
effect taken into account (dashed line, see text). The measured Fly’s Eye profile, which fits 
the profile of a nucleonic primary, is shown by the data points (Halzen and Hooper 2002). 



Fig. 9. — Two OWL satellites in low-Earth orbit observing the flourescent track of a giant 
air shower. The shaded cones illustrate the field-of-view for each satellite. 
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light in the near UV (300 - 400 nm range) that is isotropically emitted by nitrogen molecules 
that are excited by the charged shower secondaries at the rate of ~4 photons per meter per 
particle. Its estimated aperture is 1000 km 2 sr at 100 EeV after inclusion of a 10% duty cycle 
(Sokolsky 1998) . 

The southern hemisphere Auger array is expected to be on line in the near future. This 
will be a hybrid array which will consist of 1600 particle detector elements similar to those 
at Havera Park and three or four flourescence detectors. Its expected aperture will be 7000 
km 2 sr for the ground array above 10 EeV and ~ 10% of this number for the hybrid array. 

The Telescope Array will will consist of eight separate flourescence detecting telescope 
stations separated by 30 km. Its expected aperture will be 8000 km 2 sr with an assumed 10% 
duty cycle. 

The next big step will be to orbit a system of space-based detectors which will look 
down on the Earth’s atmosphere to detect the trails of nitrogen flourescence light made by 
giant extensive air showeis. The Orbiting Wide-angle Light collectors (OWL) mission is 
being proposed to study such showers from satellite-based platforms in low Earth orbit (600 
- 1200 km). OWL would observe extended air showers from space via the air fluorescence 
technique, thus determining the composition, energy, and arrival angle distributions of the 
primary particles in order to deduce their origin. Operating from space with a wide field- 
of-view instrument dramalically increases the observed target volume, and consequently the 
detected air shower event rate, in comparison to ground based experiments. The OWL 
baseline configuration will yield event rates that are more than two orders of magnitude 
larger than currently operating ground-based experiments. The estimated aperture for a 
two-satellite system is 2.5 < 10 5 km 2 sr above a few tens of EeV after assuming a 10% duty 
cycle. 

Figure 9 illustrates two OWL satellites obtaining stereoscopic views of an air shower 
produced by an ultra-high energy cosmic ray. With an approximate 10% duty factor, OWL 
will be capable of making accurate measurements of giant air shower events with high statis- 
tics. It is expected to be able to detect more than 1000 showers per year with E > 100 EeV 
(assuming an extrapolation of the cosmic ray spectrum based upon the AG AS A data). 

The European Space Agency is now studying the feasibility of placing such a light collect- 
ing detector on the International Space Station in order to develop the required technology to 
observe the flourescent trails of giant extensive air showers, to make such observations, and 
to serve as a pathfinder mission for a later free flyer. This experiment has been dubbed the 
Extreme Universe Space Observatory (EUSO) (see paper of Livio Scarsi, these proceedings, 
for more details). Owing to the orbit parameters and constraints of the International Space 
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Station, the effective aperture for EUSO will not be as large as that of a free flyer mission. 

A recent compendium of papers on observing giant air showers from space may be found 
in Krizmanic, Ormes and Streitmatter (1998). 


12. A Cloudy Present - A Bright Future 

As of this writing, there is a disagreement in the trans-GZK event rate between the 
AG ASA and HiRes experimental groups. Thus, we are uncertain about the observational 
situation. The prospect of new physics and new astrophysics at ultrahigh energies has 
produced a plethora of theoretical ideas and papers. Indeed, if there are significant numbers 
of ultrahigh energy events above 100 EeV, and especially above 300 EeV (which would rule 
out the heavy nucleus scenario) many of the theoretical models presently proposed could be 
ruled out. This situation might then call for radically new physics such as would involve 
violation of Lorentz invariance. 

New and more powerful observational techniques are called for to obtain significantly 
large numbers of giant air shower events to analyse in order to accurately determine the 
flux and energy spectrum of trans-GZK cosmic rays. The Auger ground array is starting 
operation. New space experiments, EUSO and OWL, have been propsed. This author hopes 
that they be built and flown and will provide the needed information. Such experiments have 
the potential of breaking through to new insights about the basic nature of the universe. 
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